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We report on the experimental observation of scaling in the time evolution following a sudden 
quench into the vicinity of a quantum critical point. The experimental system, a two-component 
Bose gas with coherent exchange between the constituents, allows for the necessary high level of con¬ 
trol of parameters as well as the access to time-resolved spatial correlation functions. The theoretical 
analysis reveals that when quenching the system close to the critical point, the energy introduced by 
the quench leads to a short-time evolution exhibiting crossover reminiscent of the finite-temperature 
critical properties in the system’s universality class. Observing the time evolution after a quench 
represents a paradigm shift in accessing and probing experimentally universal properties close to a 
quantum critical point and allows in a new way benchmarking of quantum many-body theory with 
experiments. 


Scaling laws and symmetries are at the foundations of 
modern science. They allow putting phenomena as dif¬ 
ferent as opalescent water under high pressure, protein 
diffusion in cell membranes [T], and early-universe infla¬ 
tionary dynamics 00 on the same structural footings. 
Typically, scaling is observed in thermal equilibrium or 
in relaxation dynamics close to equilibrium 00 while 
recently the scaling hypothesis has been extended to far- 
from-equilibrium dynamics BUMS- 

For spatially extended systems, it is natural to ask how 
its different parts are correlated with each other. Close 
to critical configurations, the essential physics is typi¬ 
cally captured by a single parameter s measuring the 
distance to criticality, and a universal function. General¬ 
izing this to non-equilibrium quantum dynamics implies 
that the universal function explicitly includes time evolu¬ 
tion. For example, given a time-dependent characteristic 
length scale £, scaling implies a self-similarity relation 

as~ uz t-,se) = s-^(t-,^ ( 1 ) 

with positive scaling factor s and critical exponents v, z. 
This relation reflects that the spatial structure for fixed 
system parameter 5 at a given time is the same as the 
structure at different e, at suitably rescaled times. 

Our experiments reveal such scaling behaviour in a ru¬ 
bidium condensate in a quasi one-dimensional configura¬ 
tion. The atoms are in two hyperfine states, which we de¬ 
note as |t) = \F = 2 ,rriF = —1) and |t) = |F = l^rnp = 
1). The collisional interaction between atoms in these 
states is tuned by use of an interspecies Feshbach reso¬ 
nance usual such that the system is in the immiscible 
regime, i.e. the two components tend to minimise their 
overlap (Fig. Bo- Specifically, the interspecies scattering 
length is chosen larger than the intraspecies scatter¬ 
ing lengths and a^. In the experiment, we chose 
a = ~ 1.23. De-mixing of such kind has al¬ 

ready been observed experimentally H3EZI, and studied 


theoretically [T8H23] . Varying the strength of linear Rabi 
coupling between the two atomic species allows tuning 
across a quantum phase transition between the immisci¬ 
ble and miscible regimes. 

Here, we study the dynamics after a sudden change, i.e. 
quench, of the linear coupling, observing the time evolu¬ 
tion of the spatially resolved density patterns of the two 
components n^(y) and n^{y) along the extended axis of 
the trap, on either side of the miscible-immiscible tran¬ 
sition (Fig. Bo- Due to the repulsive interactions be¬ 
tween the individual components, the local atomic den¬ 
sity p(y) = n^(y) + ni(y) is to a very good approxi¬ 
mation constant, such that only the density difference 
n t(v) ~ n l(v) fluctuates. Consequently, the gas can 
be considered as a homogeneously distributed, coupled 
collective spin ensemble, characterised by a continuous 
angular-momentum field J (y) (Fig. [2^). The longitudi¬ 
nal extension of ^ 200 /im allows exploring the miscible 
regime where the expected length scales are of the or¬ 
der of a few microns, as well as the immiscible side, with 
expected domain sizes of a few tens of microns. The z- 
component J z (y) = [n^(y)— n^(y)\/p(y) of the ‘Schwinger 
spin’ normalised to the total density is related to the 
density difference, and spin correlations G zz (y, y'. t) = 
(J z (y)Jz(y'))t between different points y,y' can be de¬ 
termined from the density patterns (Fig. [TJd) . 

The long-wavelength dynamics of our spin fluid is 
given, to a good approximation [24], by a translation- 
ally invariant nonlinear XXX-type Heisenberg Hamilto¬ 
nian density [25], 

u = [\d y 3\ 2 /4 + nj x - n c j*/2]p/2. ( 2 ) 

The first term represents the Heisenberg spin coupling. 
The term QJ X provides the local coupling of the two 
components, in analogy to an effective magnetic field 
acting transversely to the spin J z , with strength given 
by the Rabi frequency D of the linear coupling. In ad- 
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FIG. 1. (Color online) Details of the experimental system, (a) Phase diagram, distinguishing miscible and immiscible phases 
of the elongated degenerate Bose gas of rubidium atoms in F — 2 (blue) and F — 1 (red) hyperfine states. The state of the 
system is controlled by linear coupling of the levels, with Rabi frequency fl, and by tuning the collisional interaction between 
atoms in the hyperfine states, quantified by the relative strength a — °f inter- and intra-species scattering lengths 

(experimentally fixed to a « 1.23). A quantum phase transition occurs at Q c = pg(a — 1), with ID atom density p and 
intra-species coupling constant g. (b) The system is initially prepared far in the miscible regime, and then U quenched close 
to Q c . After different evolution times the two species are absorption imaged. Snapshots of the patterns emerging on either 
side of the transition are shown, with corresponding normalised density imbalance (n^(y) — n\,(y))/p and density correlation 
functions between spatially separated points y and y' . The correlations on the miscible side exhibit decay on a characteristic 
length scale, while oscillations on the immiscible side reflect domain formation as seen in the density. 


dition, a ‘single-ion anisotropy’ [26; term J% appears 
which results from the local collisional interaction be¬ 
tween the two components, with Q c = pg(a — 1) propor¬ 
tional to the tunable interaction strength between the 
spins. As the scattering lenghts of the respective Rubid¬ 
ium hyperfine scattering channels are very close, we take 
= g/(fud_ l), such that a = huo±a^/g. This 
system reveals, at zero temperature, a quantum phase 
transition at U = Q c where the Rabi-induced mixing of 
the components cancels the effect of the interspecies scat¬ 
tering, with the order parameter given by the magneti¬ 
zation (J z ). Specifically, for a strong effective magnetic 
field U > Q c the spins in the ground state are polarised 
in the x direction, i.e. the two components can not spa¬ 
tially separate although the bare system (ft = 0) is phase- 
separating. This is confirmed experimentally (Fig. ER 
A generic scaling hypothesis which includes dynamics 
out of equilibrium implies that the spin-spin correlations 
(■ Jz(y)Jz{y'))t,e = G zz {y - y',t;e) after a sudden quench 
of the linear coupling obey 

G zz {s~ v y, s-" z t- se) = s~^G zz (y, t; e), (3) 

where 5 = Q/Q c — 1 is determined by the final effec¬ 
tive magnetic field U after the quench, v and z are 
critical exponents, and g is known as anomalous expo¬ 
nent. In the experiment we reach this out-of-equilibrium 
regime by initially preparing the system with a fast 7r/2 
microwave-radiofrequency pulse in a J x spin state which 


is the ground state of the system in the infinite-linear- 
coupling limit (0 Q c ). Then, the intensity of the 
radiofrequency field is quickly reduced, switching 5 to its 
final value. Adjusting the linear coupling during the fol¬ 
lowing evolution compensates for the change of Q c due to 
the loss of particles which was independently determined. 

The correlations G zz {y,t]e) developing on the misci¬ 
ble side (e > 0) are shown in Fig. [2 Jd, in comparison with 
homogeneous Bogoliubov-de Gennes theory predictions 
[27], averaged over the density inhomogeneity in the trap 
and convoluted with the optical point spread function 
of the imaging system. Fitting an exponential to the 
short-distance fall-off of the observed correlation func¬ 
tions we extract a correlation length £(t;e) which shows 
near-linear growth after the initial quench (Fig. ED- The 
Bogoliubov prediction (solid lines) qualitatively repro¬ 
duces this rise as well as the oscillations seen for larger 
e. The damping of the oscillations seen at smaller e is 
attributed to effects of the transverse trapping potential. 

We extract the maximum correlation length for differ¬ 
ent e within the first 12 ms after the quench, where this 
observable is still weakly affected by the atom loss. Using 
the theoretically expected exponent of 0.5 for rescaling 
the correlation functions at a fixed time (t = 12 ms) they 
all fall on a universal curve (Fig. upper panels). Ex¬ 
tracting characteristic length scales as indicated, we find 
scaling according to Eq. 0 (lower panels of Fig. [§i). The 
exponent extracted from a linear fit of £ on a double-log 
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FIG. 2. (Color online) Time evolution of correlations after 
quench to the miscible side of the quantum phase transition, 
(a) Two-component gas as a coupled collective spin ensem¬ 
ble. The spatially resolved density difference between F — 1 
and F — 2 allows the extraction of the local ^-component of 
the collective spin vector J (y). (b) Measured spin correlation 
function G zz (y,y',t ) = {J z (y)Jz(y'))t at three different times 
after a quench to £ = Q/Q c — 1 = 0.17, showing build-up 
of spatial correlations. Solid black lines show Bogoliubov-de 
Gennes mean-field predictions, (c) Time-evolution of the cor¬ 
relation length £(£; e), for three different e, deduced by fitting 
an exponential to the short-distance fall-off of the extracted 
correlation data. The Bogoliubov evolution (solid lines) re¬ 
covers the initial near-linear rise of £, with slope given by the 
speed of (spin-wave) sound, and captures the maximum cor¬ 
relation length at a characteristic time depending on e. The 
predicted oscillatory behaviour is experimentally observed for 
larger e, while the maximum correlation length reached at 
short times is robustly detected in all cases. Dashed lines serve 
as a guide to the eye, marking the predicted first maxima. 
For the scaling analysis of the data, the correlation lengths 
are compared at a fixed time, see Fig. [3] 

scale is v = 0.51 ± 0.06 (Fig. UK lower right panel). The 
scaling exponent is robust with respect to varying the 
range of the exponential fit of the correlation functions 
and the way of extracting a characteristic length after 
the initial linear rise. In the immiscible regime, e < 0, 
we choose the domain size L& and find an exponent of 
v = 0.51 =t 0.04 (Fig. If, lower left panel). As a result, 
in both, the miscible and immiscible regimes, we find self 
similarity under rescaling y -F e u y with v = 1/2, with 
correlations following different universal functions. 

To obtain the dynamical exponent z we analyse the 
observed time dependence of the correlation functions, 
focusing on the point of time r where £(t) levels off after 
the initial linear rise (Fig. [3 Jd, upper panel). For this we 
determine the crossing point of two independent linear 
fits to the short and long-time behaviour. The resulting 


dependence of r(e) is compared, on a double-log scale to a 
power law with exponent vz = 1/2 (Fig. If, lower panel, 
dashed line). The deviation of our data from this power 
law for large e can be understood within Bogoliubov 
theory which predicts £(£;£) = (2raD c £) _1 / 2 | sin(A(e)t)| 
m- The time ti j£ of the first maximum (solid line in 
Fig. |f, lower panel), related to the characteristic time r, 
is inversely proportional to the gap, i.e., zero-momentum¬ 
mode frequency A( 5 ) = Q c y/e(e + 1). 

On either side of £ = 0, the critical exponents extracted 
from our data are consistent with Bogoliubov-de Gennes 
mean-field predictions. To reveal the limitations of mean- 
field theory and study the effects due to the excitation 
of total-density fluctuations, we have performed semi- 
classical simulations of coupled Gross-Pitaevskii equa¬ 
tions for the two components (Fig. [4|. We define a cor¬ 
relation length £0 in terms of the zero-momentum limit 
of the Fourier transform of the spin-spin autocorrelation 
function, G zz (k : t;e) = / dyexp{—iky}G zz (y, t\ e) [123] . 
For £ > 0.1, we find very good agreement of the computed 
correlation length at the first maximum (blue points) 
with mean-field scaling (solid black line). The numerical 
data shows that, for £ < 0.05, the extracted correlation 
length scale at the first maximum deviates from a sim¬ 
ple mean-field power law and saturates to a finite value 
at vanishing e. This £-dependence of £0 shows universal 
cross-over behaviour reminiscent of an equilibrium one¬ 
dimensional (ID) Ising system [28]. The spontaneous 
breaking of the discrete Z 2 symmetry under J z -F —J z 
allows for spin-t/i domain formation and thus, according 
to the standard theory of critical phenomena, the transi¬ 
tion in our system belongs to the Ising universality class 
m • The analysis of the numerically observed crossover 
behavior requires a discussion of non-perturbative correc¬ 
tions to the experimentally observed mean-field scaling. 
As seen in Fig. [4] our current experiment touches on the 
regime where non-perturbative corrections start to be¬ 
come important. 

In conclusion, we have demonstrated the possibility of 
probing universal properties close to a quantum critical 
point by quenching the system out of equilibrium and 
observing the short-time evolution long before equilibra¬ 
tion. With that, our experiment opens a new path to 
study universal properties building on phase coherence in 
closed many-body systems. This is essential for bench¬ 
marking quantum many-body theory, moving towards a 
quantum simulator of universal critical phenomena. 

Acknowledgements. The project was initiated in ex¬ 
tensive discussions with Jacopo Sabbatini and Matthew 
Davis in the context of the Kibble-Zurek mechanism in 
ID spin systems. We acknowledge the important in¬ 
put from Isabelle Bouchoule on the Bogoliubov descrip¬ 
tion of the observations, and thank Tilman Zibold for 
support on the experiment. We thank J. Berges, S. 
Erne, F. Essler, and J. M. Pawlowski for discussions. 
This work was supported by Deutsche Forschungsge- 






































4 


(a) space rescaling 




(b) time rescaling 




FIG. 3. (Color online) Spatial and temporal scaling of the spin-spin correlations, (a) Spatial correlations after quenches to 
different proximities e from the critical point (color-coded), on the immiscible (left panels, at t — 39ms) and miscible side 
(right, t — 12 ms). Top row: The correlation functions, under a rescaling y —y ye 1 ' of the distance dependence with mean-held 
exponent v — 1/2, fall on a universal curve. Bottom row: ^-dependence (double-log scale) of the characteristic length scales 
deduced from the correlation functions. The straight lines reveal values for the critical exponent v — 0.51(4) on the immiscible 
and v — 0.51(6) on the miscible side of the transition, (b) Temporal scaling of the spin-spin correlations. The characteristic 
time t for different £ is obtained as the intersection point of the linear fits to the initial rise of £(£) (grey symbols for s = 0.1) 
and to the behaviour after £ deviates from this rise. The procedure of determining the intersection is exemplarily shown in the 
upper panel for e = 0.23. In the lower panel we compare the extracted r(e) to a mean-field scaling with vz — 1/2 (dashed line) 
and to the Bogoliubov prediction r ~ 1/A, with gap A(e) = Pl c \Js(e +1), also applicable at larger e > 1 (solid line). 



FIG. 4. (Color online) Scaling analysis at short times after 
the quench. Correlation length £o(£; e) on the miscible side 
of the transition, at the time of the first maximum of £o af¬ 
ter the quench. The solid black line marks the Bogoliubov 
prediction £Bog = \fh/ (2m£l c e) also shown in Fig. [ 3 ^ 1 . Closer 
to the critical point, the open coloured symbols show results 
of semi-classical simulations of the quench dynamics. In the 
experimental range (e > 0.1) simulation data and Bogoliubov 
mean-field prediction agree. For e < 0.1, the simulations 
demonstrate a deviation from the mean-field power law and 
saturation for £ —0. 
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SUPPLEMENTARY MATERIALS 


|2j, for = an = a = g/2huj±_ (in units where h= 1), 


Spin-fluid representation of the two-component 
Bose gas 


Model. The microscopic Hamiltonian for two interact¬ 
ing Bose fields (h/, j G {t, i}, in one spatial dimension 
(ID), with a linear Rabi coupling Q and a zero-point de¬ 
tuning S is H = Hq + H C p\ + -Hint 5 with the quadratic and 
quartic terms 


H 0 = 



X ( d 2 

2 m v 


V(y) 


*3, 


Hc P i = ^Jdy + h.c.) + - $}$+)], 

]T 9ij [dyQiQ&fr. (4) 

1 i,3= td J 


Here, m is the atomic mass, V(y) = muj 2 y 2 /2 the lon¬ 
gitudinal trapping potential, and < 7 ^ = 2huo_\_aij are the 
effective ID couplings, with intra-species 3D scattering 
lengths Off, an an d inter-species scattering length oqq. 
In the experiment, a = &tt/-\/ a tt a U is tuned to a = 1.23 
by means of the inter-species Feshbach resonance at 
9.10 G. For Q = 0, S = 0, and equal intra-species cou¬ 
plings = an , a non-vanishing cross-coupling, a ^ 1 , 
leads to a deviation from a fully U (1) x SU (2)-symmetric 
theory. Adding the linear Rabi coupling, 9^0, leaves 
a residual symmetry, which is broken if S ^ 0 or 

Spin-fluid representation. The model © of coupled 
Bose gases can be interpreted to describe one gas of 
atoms in two spin states. We term this system a co¬ 
herent spin fluid, as each volume element of it carries 
angular momentum proportional to the local density p 
of atoms, generalizing the concept of classical relativis¬ 
tic spin fluids pQ. This fluid is not to be mixed up 
with the disordered state of a spin liquid, occurring, e.g., 
in lattices giving rise to frustrated moments. In the 
case that mode occupation numbers are large as com¬ 
pared to 1, the Bose field operators can be replaced by 
classical complex fields, <f >j —)> (f)j. The resulting clas¬ 
sical Hamiltonian is separated into internal degrees of 
freedom (DoF) and hydrodynamic ones by use of nor¬ 
malised Schwinger angular momenta J a = p _1< F-cr?<Iq 
with a £ {x, y , zj (sum over repeated indices is im¬ 
plied), with Pauli matrices a a and p = p^ + Hence, 
Jx = Zp -1 i/PiPi cos ip, Jy = 2p~ 1 ^/p^plsiiup, J z = 
p~ l {p\ — Pi) , normalised to |J| = 1. Here, <p = <p-\ — <pi is 
the relative phase of the fields <pj = ^flpjexp (upj). The 
3-projection J z of the spin density describes the relative 
density difference of the fluids. In terms of spin densities 
and total DoF, p and 0 = ipi + tpi, the Hamiltonian reads 


H =\ Jdy^rn 1 (d yy /p) 2 + £^\d y d\ 2 + p(V + mv 2 s ) 

2 a 

+ 9P 2 + ^fl( a — 1) [l — (Jz) 2 ] + FtpJ x + SpJz\ • (5) 


Here, the effective velocity of the fluid, v e s = 
pl y l)/p, defined in terms of the velocities = 
of the spin components, is defined as [2] 


^eff = 


2m 


dy e 


Jz(JydyJx J X dy jy) 


[1 - (Jz) 2 } 


(pt v t + 

dzVi/m 


(6) 


The Hamiltonian © expresses the energy of the two- 
component Bose gas as that of a ID spin field J (y) carried 
by a fluid with density p and conserved current j = pv e s. 
For a fluid at rest, i.e. p = const, and v G s = 0, the spin 
system assumes the form of a nonlinear sigma model [2] 
or XXX-type Heisenberg chain with single-ion anisotropy 
~ J z in a transverse and/or longitudinal magnetic field 
[3], with an additional chiral-field term oc v 2 H . This, as 
well as the anisotropic terms oc [1 — ( J z ) 2 ] and oc J x ,Jz 
break the 0(3)-symmetry of the spin system. 


Quantum phase transition 

Transition between miscible and immiscible phases. To 
set the stage, before we describe the non-equilibrium 
dynamics realized in the experiment, we briefly sum¬ 
marize the equilibrium phase structure of the system. 
The non-linear Heisenberg model (J5| for the spin DoF 
(p = const., t’eff = 0) represents a classical isotropic 
Heisenberg ferromagnet. The homogeneous (V = 0) Bose 
gas © of two linearly uncoupled (Q, 5 = 0) components 
possesses two distinct ground states depending on the 
choice a mg. In the immiscible regime, a > 1, the 
term oc (a — 1)[1 — ( J z ) 2 ] gives a positive contribution 
to the total energy for all possible spin configurations, 
except when J z (y) = d=l. Spontaneous breaking of the 
remaining discrete Z 2 symmetry under J z —> — J z leads, 
in a non-equilibrium or thermal system, to the formation 
of domains in the spin density J^, i.e. oppositely signed 
patches of J z with J z ~ ±1, similar to the ferromag¬ 
netic classical Heisenberg model [?, [8]. In between the 
domains, due to the chiral term cx v 2 s , the spin varies in 
a way that rotations around the J z axis are suppressed. 
On the contrary, in the miscible regime, a < 1, spin 
configurations with J z (y) = 0 are preferred energetically, 
and thus in the zero-temperature ground state one has 

(H(y) = Pi(y)- 

Modulational instability. At zero temperature and for 
Q = S = 0, the system is immiscible for the experimen¬ 
tal choice a > 1. It can be rendered miscible again by 
introducing an additional linear Rabi coupling Q, with 
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Cl > Cl c = pg(a — 1) = — pg s . A Bogoliubov anal¬ 
ysis i of the homogeneously mixed system gives, for 
= a U = 0? two decoupled excitation branches 

w+(«, a) = Cl c ^/(a + l)(<a - l) _1 ft 2 + ft 4 , (7) 

cj_(ft,£) = 5(5 —|— 1) —|— (2e l)ft 2 -j- ft 4 (8) 

of density and spin waves, respectively, where ft 2 = 
k 2 /[2mCl c \, and 6 = (Cl — Cl c )/Cl c is the relative prox¬ 
imity to the critical coupling Cl c . For unequal couplings, 

| an ~ a ttl ^ a hb the two branches hybridise, develop¬ 
ing avoided crossings [9|. For a > 1 and — 1 < 5 < 0, 
the spin-wave frequency c is imaginary for modes with 
momenta ft < k c = \f^e. These modes hence be¬ 
come unstable, resulting in spin waves with growing am¬ 
plitude P and leading to a spatial separation of the 
two components. The most unstable mode is that with 
ft 2 = ft 2 = max{-e — 1/2,0}. In the presence of the 
full non-linearity of the purely repulsive two-component 
Gross-Pitaevskii model this classical modulational insta¬ 
bility implies small initial density fluctuations on top of 
the homogeneous mixed system to grow into fully po¬ 
larised areas delimited by sharp boundaries with a thick¬ 
ness given by the spin healing length = y/2/k c . 

In the miscible regime, e > 0, the spin-wave modes 
acquire a gap cj_(fc = 0) = — Cl c \ = A(e) = 

Cl c ^/e(e + 1) in their excitation spectrum, such that the 
mixed state remains stable mm- The disappearance 
of the gap at 6 = 0 defines the critical line of a second- 
order (quantum) phase transition in the plane (Main 
text, Fig. 00- Note that the linear coupling oc Cl in the 
Hamiltonian breaks the 0(2)-symmetry in the xy- 
projection of the spin configuration space. As a conse¬ 
quence the conservation of the total polarisation N\ — N 2 
is lost. 

Experimental procedure 

We prepare a condensate of N ~ 4 x 10 4 87 Rb atoms 
in the |}) = \F = 1 ,nriF = 1) hyperfine state of the 5 Si/ 2 
manifold, in an elongated optical trapping potential with 
(cj y ,(jj±) = 27r x (1.9,128) Hz, resulting in a peak density 
nip — 230/im -1 . The Feshbach resonance between |}) 
and |t) = |F = 2,mf = —1) at 9.10 G is used to tune 
the inter-species scattering length to ~ 120aB O hr at 
9.08 G while the intra-species scattering lengths are fixed, 
(a^,an) = (95,100) aBohr- At this magnetic field, the 
bare system is immiscible. By applying resonant two- 
photon linear coupling, the system can be rendered mis¬ 
cible. The phase transition occurs at the critical Rabi 
frequency Cl c . In mean-field and uniform-gas approxi¬ 
mation, it is given by Cl c = — pg s « 27t x 70 Hz where 
g s = fud±(a^ + aii — 2 a^) is the ID coupling of the spin 
degrees of freedom, and p the ID density. 


After the initial preparation of all atoms in |}), a 7 t/ 2 
pulse of combined microwave and radio frequency (RF) 
magnetic fields with a Rabi frequency Cl « 27t x 340 Hz, 
creates an equal superposition of | }) and 11) • Subse¬ 
quently, the phase of the radio frequency is switched by 
7 t/2 (37t/ 2 for e < —1), and Cl is quenched to a final 
value above/below the critical value Cl c . After the follow¬ 
ing evolution time, the two components are sequentially 
detected using in-situ absorption imaging at 9.08 G. 

Due to the proximity of the Feshbach resonance, atoms 
are lost with a 1/e lifetime of ~ 30 ms, and the density 
decreases with time. The resulting change of Cl c was com¬ 
pensated by dynamically adjusting the Rabi frequency Cl 
during the evolution time. All given Cl refer to the initial 
value. Other effects of atom loss such as the change in 
the spin healing length remain. To ensure resonance of 
the two-photon coupling, the independently determined 
AC Zeeman shift resulting from the detuning of 200 kHz 
to the intermediate |F = 2 ,mp = 0) state is compen¬ 
sated by adjusting the frequency of the RF. The average 
of the density dependent mean-field shift oc (g^ — g^)p 
is compensated in the same way. 

The longitudinal extension of ~ 200 pm. allows explor¬ 
ing the miscible regime where the length scales are of the 
order of a few microns as well as the length scales of a 
few tens of microns on the immiscible side. The trans¬ 
verse extension of the atomic cloud is ~ 2 pm (chemical 
potential p& 2n x 300 Hz) and is comparable to the spin 
healing length £ s (Cl) = h/^m(g s p + hCl) = y/h/mCl c s 
for |^| ^ 1. Thus the system close to quantum critical¬ 
ity is effectively one-dimensional for the spin degree of 
freedom. 

In the analysis of the images we select the center of 
the cloud to reduce effects of cloud inhomogeneity and 
position. While the pixel resolution of the CCD cam¬ 
era is 0.41 /im, the resolution of the imaging optics is 
larger by about a factor of three. Hence, we group pixels 
by three and calculate correlation functions on a spatial 
grid of 1.23 pm. The experimental data is analysed with 
respect to the linear density profiles p\^(y) of each com¬ 
ponent from which the spin profile J z (y) and its longitu¬ 
dinal correlations G zz (y,y' ,£) = (J z (y)Jz(y'))t between 
different points //, y' along the axis of the trap are com¬ 
puted, and averaged over ^ 20 experimental realisations. 
We extract the correlation length £(t;e) by fitting an ex¬ 
ponential exp{— y/[y/2£(t; e)]} to the small-7/ fall-off of 
the correlation function G zz (y, 0, £; e) normalized to 1 at 
y = 0. In the data analysis, images with an integrated 
population imbalance | — N^\/[N^ + IV}] > 0.2 are not 
taken into account, as well as images with large imaging 
noise. To reduce the effects of the inhomogeneity we sub¬ 
tract the spin profile, smoothed by a Gaussian filter with 
a width of 80 /im, from the measured profile. The width 
ensures that structures and fluctuations on the scale of 
the correlations to be measured are conserved while the 
profile gets centered around zero imbalance. This pro- 








cedure corresponds to applying a high-pass filter to the 
spin profile. 

Critical dynamics following a quench into the 
vicinity of the phase transition 

Critical scaling in equilibrium. Before we continue with 
the theory of the dynamics induced in the experiment, we 
briefly recapitulate equilibrium mean-field scaling prop¬ 
erties. We continue to consider the case a > 1. In 
the vicinity of the quantum phase transition at Q c = 
pg(a — 1), correlations show scaling behaviour with re¬ 
spect to the relative distance 5 to the critical point. Equi¬ 
librium renormalisation-group theory predicts the scaling 
of the correlation length near criticality in terms of the 
critical exponent v as defined in Eq. 0 - The mean-field 
prediction for v can be derived within semi-classical Lan¬ 
dau theory as follows. Neglecting fluctuations of the total 
density p and phase 0, the Hamiltonian 0 can be writ¬ 
ten in terms of Euler angle representation of the spin, 
J x = cos#cos(^, J y = cos 9 sin<p, and J z = sin#, as 

H = l f dy {f^ + ( d v°) 2 } 

+ 2f2cos# cos(p + f2 c cos 2 (9) 

with —7r/2 < 9 < 7r/2, — it < cp < 7r, and ip and 9 being 
mutually dual variables while d y (\) = Uq = d t 0 is the 
canonical conjugate of 0. For Q c = 0, this model bears 
an 0(2) symmetry around the J x axis and undergoes 
a quantum phase transition in the Kosterlitz-Thouless 
class El. Choosing Q c > 0, the energy density exhibits 
the possibility of spontaneous Z 2 symmetry breaking at 
5 = 0 which becomes apparent by expanding around 6 = 

0, = 7T, 

-(1 + 2 e)+e8 2 + + 0(d 6 ) . 

(10) 

While for 5 > 0 there is only one minimum at 6 = 0, the 
ground state can assume, for —2 < 5 < 0, different values 
#o = arccos(l + e). Hence, while the potential remains 
symmetric under 0 —>> —0, the ground state is two-fold 
degenerate with respect to the spin orientations J z = 
±^/\e\{2 — \e\) = zb\/^c — and the Z 2 symmetry 

is spontaneously broken. (9) or, equivalently, (J z ), is the 
order parameter of the transition. 

The inverse coherence length squared is proportional 
to the ‘mass’ parameter in the potential, i.e. the second 
derivative of H po t with respect to 9 at the value of the 
order parameter. In mean-field approximation, it follows 


from the potential ( fl()| ) that 

2 / £ (5 > 0, symmetric phase), 

^ 0(1 { \e\(2 — |e|) (—2 < 5 < 0, broken phase). 

( 11 ) 

Together with Eq. 0, this shows that, in both phases, 
v = 1/2. To mean-field order this result is independent 
of the dimensionality D m d + 1. 

Dynamical exponent Out of equilibrium, universality 
classes are defined, in addition, by the dynamical expo¬ 
nent z [12]. From Eq. (J8| we infer, for 5 <C 1 and k <C 1, 
that 

u)-(sk, se) = s z uj-(k, 5 ). (12) 

and therefore z = 1 at mean-field level. Our numerical 
study of the system’s response to the experimentally ap¬ 
plied quench, discussed in the following, shows that the 
above mean-field scaling is approximately valid in the 
non-equilibrium situation of the experiment. 

Initial state and quench. In the experiment, the two- 
component gas is prepared in a quasi-condensate, with 
interactions tuned to a fixed a 1.23 > 1, equal occu¬ 
pation of the modes j — t and and zero relative phase, 
i.e. fully polarised along the positive spin x-axis. In this 
state, the relative DoF can, to a good approximation, 
be considered in the zero-spin-temperature ground state 
while the centre DoF are equilibrated at a total temper¬ 
ature on the order of T ~ 30 nK. The system is then 
subject to a quench in D close to the value where the 
above described quantum phase transition occurs in the 
equilibrium setting. 

Bogoliubov description of the time evolving spin cor¬ 
relations. We briefly summarize the results of a Bo- 
goliubov mean-field analysis of the spin degrees of free¬ 
dom, for further details cf. Ref. [2- For convenience we 
introduce the dimensionless position y = ?/(2mO c ) 1 / 2 , 
time t = £D C , linear density p = p/(2mD c ) 1 / 2 , mo¬ 
mentum k = /c/(2m^ c ) 1 / 2 , mode frequency uj-(k,s) = 
cj_(ft,£)/f2 c , gap A = A/D c = ^ 5(5 + 1), spin-wave 
speed of sound c s = (2m/D c ) 1 ^ 2 c s = %/2e + 1, and tem¬ 
perature T = kBT/D c . 

The Fourier mode expansion of the angle fields is 
<p(y,t) = (2\/p) _1 X)« ; [/^(y)^ ex P{-^«i} + h -c-L and 
9(y,t) = (2^/p)- 1 J2 K [U^(y)b K exp{-ioj K t} + h.c.\. For 
a homogeneous system one finds the Bogoliubov-de 
Gennes mode functions f±(y) = Z -1 / 2 [o)_(ft, £)/(k 2 + 
5 + l)] ±1 / 2 exp {iny}. L is the spatial length, and the 
quasiparticle operators obey [6^,6^,] = S KK /, [b K ,b K r] = 
0. The experimental range e G {0.1...1} implies 
A G {0.33... 1.4}, K g G {0.3... 0.8}, c s G {1.1... 1.7} 
and thus a relatively small near-linear sound-wave 
regime. The range of c s implies sound speeds c s G 
{0.44 ... 0.68} pm/ms. 

In the experiment, the spin, pointing initially into 
the negative z-direction, is rotated into the x-direction, 
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i.e. to (0,(p) = (0,0). After this, the Rabi coupling Q is 
quenched close to the critical coupling, 6 ~ {0.1... 1}. 
As a consequence of the initial tt/2 rotation of the spin, 
the initial spin and relative phase fluctuations can be 
considered to be Gaussian and delta-correlated in space: 
G Z z{y,y'-,t = 0) = {J x {y)Jx{y '))o - {9{y)9{y')) 0 = 
(2p)~ 1 5(y - y') = (<p(y)<p(y')) o- 
The local initial fluctuations of 0 and ip can be used to 
compute the time evolution of the spin-spin correlation 
functions, 


968 (y,t ) = (9(y,i)9(0,t)) 
= (2/3) -1 S(y) + j 


,%Ky 


sin 2 (o)_(/^, e)i) 

K 2 + £ 


(13) 


Using the approximation u)_ ~ A + c s k + fir in Eq. (13) 


one can perform the momentum integral analytically. 
The result gives an exponential decay with an oscillatory 
behaviour added, which dominates the function at early 
times and short distances (Main text, Fig. [ 2 ^ 1 ), goo(y , t) = 
(2p)~ 1 [S(y) + (4 v / 5) _1 (exp{- A /5|^|}+oscillatory part)]. 
We note that the extension of the cloud in the trap with 
longitudinal frequency uj z ~ (2i r) 2 Hz allows an estimate 
of the relative size of the exponentially decaying part 
to the delta-noise term to ^uj z /Q c £/4: = {0.04... 0.13}. 
Hence, the correlation function in Gaussian approxima¬ 
tion is 


G zz (y,t) = (, J z (y)J z (0)) t 

= sinh ([g$e(y, t) + gee(~y,t)}/2 ). 


(14) 


In the experiment, p ~ 210, such that, away from y = 0, 
G zz {y, t) is well approximated by [go$(y, t)+gee{-y, t)]/2 . 

To compare our experimental data with Bogoliubov 
theory as shown in Fig. [2 Jd and c, we evaluate the full 
correlation function G zz (y,t ) during the time evolution 
after a quench from an initial value Hq = 10 U c to differ¬ 
ent final Q close to the critical value, taking into account 
a 1/e life time of the subsequently decaying system of 
42 ms. We bin the resulting function corresponding to 
the CCD resolution of 0.41 g m and average over three 
neighboring bins with equal weight to take into account 
the finite resolution of the imaging system. The corre¬ 
lation length is extracted by an exponential fit to the 
normalized correlation function, as done for the experi¬ 
mental data. 

Numerical scaling analysis. The correlations discussed 
above are universal to the extent that, neglecting the os¬ 
cillatory part, they are characterised by a single relevant 
parameter. In our numerical analysis we take this pa¬ 
rameter to be the correlation length £ n which is defined, 
after subtraction of the S(y) contribution of the initial 
state, as the square root of the integral of the normalised 
correlation function 2pgoo(y,t;e) over all y. According 
to Eq. (13) the mean-field value of £ n is 


Ibo S (t;e) = e 1/2 | sin(At)| 


This expression for the oscillating correlation length re¬ 
produces qualitatively the behaviour seen in the data and 
represents the mean-field universal scaling function obey¬ 
ing the scaling relation £,Bog{ s ~ uz t] s ^) = s^^BogUj^) 
with v — 1/2, z = 1. 

Since the critical behaviour we are interested in is 
dominated by the large infrared (IR) mode occupation 
numbers, we can make use, for studying the non-linear 
dynamics of our system beyond the mean-field approxi¬ 
mation, of semiclassical field simulations, also known as 
the Truncated Wigner approach [13 HU, which is non- 
perturbative and goes far beyond the mean-field treat¬ 
ment discussed above. As a result, the dynamics of these 
modes can be represented by an ensemble of field tra¬ 
jectories which are propagated according to the classical 
equations of motion derived from Eq. Q, 

i dt4>^ = [Ho + g( |</q| 2 + <a|</q| 2 )] , (16a) 

i dt4>i = [Ho + g( |</q| 2 + <^|0tl 2 )] 5 (16b) 


where Ho = —dy /(2 m) + V(y). Initial field configurations 
0td(^, ^o) are sampled from a Gaussian Wigner distribu¬ 
tion, which is positive everywhere and takes into account 
initial quantum fluctuations. At a given evolution time, 
correlation functions are obtained by averaging the corre¬ 
sponding observable over the ensemble of sampled trajec¬ 
tories. For our simulations we choose V = 0. As before, 
our simulations are done in a fully one-dimensional ge¬ 
ometry. 

In accordance with the experimental procedure, we 
consider an ensemble of initial states at zero spin tem¬ 
perature, showing ground-state Gaussian fluctuations 
around a fully J^-polarised configuration. Other than in 
experiment we assume the inter-species coupling to van¬ 
ish initially, a = 0. In our simulations, the gap parameter 
A is suddenly quenched, at time t = 0, from a large pos¬ 
itive value to positive values close to zero, while simulta¬ 
neously switching the inter-species coupling from zero to 
cy = 1.23. We then measure the time-dependent spin-spin 
correlation function G zz (y,t]£) = (Jz(y)Jz{ty)t,e as a 
function of y and infer, in accordance with the above def¬ 
inition for the mean-field case, a spin correlation length 
£ n (£, e) from the zero-momentum value of the correspond¬ 
ing Fourier spectra. 

Our numerical evaluation shows, in quantitative ac¬ 
cordance with the mean-field result (15), an initial rise 
in time and oscillations around a mean value. At long 
times, the oscillations are damped, however, and the cor¬ 
relation length tends to an asymptotic stationary non¬ 
zero limit. To compare with the mean-field critical be¬ 
haviour discussed before, the time evolution of the cor¬ 
relation length £ n (£;e), evaluated at the time t = ti :£ of 
the first maximum, is studied with respect to its scaling 
in 5. Agreement with mean-field scaling is found within 
the experimentally accessed regime of 5 as discussed in 
the main text, cf. also Fig. [4j 


( 15 ) 
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